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Abstract: We show that a longitudinal gauge degree of freedom for a vector field is 
equivalent to a Pais-Uhlenbeck scalar field. With the help of this equivalence, we can 
determine natural interactions of this field with scalars and fermions. Since the theory has 
a global £7(1) symmetry, we have the usual conserved current of the charged fields, thanks 
to which the dynamics of the scalar field is not modified by the interactions. We use this 
fact to consistently quantize the theory even in the presence of interactions. We argue that 
such a degree of freedom can only be excited by gravitational effects like the inflationary 
era of the early universe and may play the role of dark energy in the form of an effective 
cosmological constant whose value is linked to the inflation scale. 
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1. Introduction 

The standard model of cosmology provides a very successful description of our Universe Q . 
However, it is based on the idea that about 96% of the energy density in the Universe stems 
from particles and fields which are not part of the standard model of particle physics and 
have never been observed in colliders. Dark matter, which amounts to about 26% of the 
Universe' content, at least has the properties of typical non-relativistic particles, but dark 
energy, the component which makes up about 70% of the energy density of the Universe 
must be endowed with very exotic properties, like a strong negative pressure, in order to 
explain the observed current acceleration of the expansion of the Universe. This most 
unexpected discovery has been awarded the 2011 Nobel prize in physics Q. 

Even if a cosmological constant has the right properties and is in agreement with 
cosmological observations, this solution is unsatisfactory from a purely theoretical point 
of view, as its associated scale (pa) 1 ^ — 10~ 3 eV is so much smaller than the natural 
scale of gravity given by M p ~ 10 18 GeV. On the other hand, one could expect this 
scale to be related to some cutoff scale coming from particle physics, which should be, at 
least, the scale of supersymmetry -E S usy > ITeV, again much larger than the value of the 
cosmological constant inferred from observations (see |J for an extensive discussion about 
the cosmological constant problem). 

This fine tuning issue of the cosmological constant has led the community to search for 
different solutions to the problem of accelerated expansion. Different dark energy models 
like e.g. quintessence or large scale modifications of gravity have been explored 

Researchers have also looked into theories with higher derivatives in the Lagrangian 
density, like e.g. the so-called f(R) theories || or the galileon field Q. Theories with 
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higher derivatives look dangerous at first sight because they typically lead to the presence of 
ghosts (particles with negative kinetic energy) or other instabilities such as the Ostrogradski 
instability |7|, ^] . This instability appears because these theories lead to higher than second 
order equations of motion, so that new degrees of freedom appear and they are usually 
ghost-like. It has been shown recently ||, that trying to eliminate the ghost by introducing 
additional constraints does not work in general because the last surviving degree of freedom 
is always the ghost. Even though without the presence of any interactions, the energy of 
a given field remains constant and the instability cannot develop, we expect that the 
coupling of such a field to other degrees of freedom or quantum effects will potentially 
provoke uncurable instabilities. 

However, exceptions to this generic rule exist in degenerate theories like the aforemen- 
tioned f(R) theories of gravity and, more recently, the galileon fields. Interestingly, the 
dark energy problem, pushes us to re-think our boundaries of what we can accept as sensi- 
ble (effective) physical theories [1C]. This is the youngest example of the cross fertilization 



of cosmological observations and fundamental theoretical physics. 

The present paper inscribes in this framework. Here we show that the action for a 
vector field having only a residual gauge symmetry, which has already been proposed as 
a candidate for dark energy |Tl]| , can be rewritten as an ordinary U(i) gauge vector field 
theory together with a (degenerate) Pais-Uhlenbeck (PU) scalar field, the field version of 
the PU oscillator first discussed in fll2|| . Even though this is a higher derivative oscillator, it 
has been shown |l3| that this model can be quantized consistently because its Hamiltonian 
corresponds to a special class of Hamiltonians that, although being non-hermitian, exhibit 
a VT symmetry that allows to construct a consistent quantum theory without negative 
energy states. Here we derives an alternative way to consistently quantize the degenerate 
PU model (at the field level) that makes use of the presence of another symmetry in which 
the PU field can be shifted by an arbitrary harmonic function. We can use this symmetry 
to restrict the physical Hilbert space or to fix the gauge and quantize only the healthy 
physical mode. Then, we include interactions with charged scalars or fermions, in a way 
which is motivated by the interpretation of this field; and we show that our restriction of 
the physical states remains intact. 

The remainder of this paper is organized as follows: In the next section we show the 
relation of the vector field action with the residual gauge symmetry and the PU field. In 
Section ^ we include interactions. In Section |3| we discuss the quantization of the theory 
and show how it could play the role of dark energy or of the inflaton. (Even if in its 
present form the model does not propose a mechanism to end inflaton.) In Section || we 
conclude and discuss further investigations which can be performed to compare this model 
with standard ACDM. 

2. The Stiickelberg trick and the Pais Uhlenbeck field 

2.1 From the Pais Uhlenbeck oscillator to the Pais Uhlenbeck field 

Let us first briefly introduce the PU oscillator \12\ and show how it arises as the Fourier 
modes of a certain higher order field theory. The PU oscillator is described by the following 
action: 

S PU = l J dt [z 2 - (to 2 + u 2 )z 2 + lo 2 u 2 z 2 ] , (2.1) 
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where 7 is an arbitrary parameter. Since the action depends on the second time derivative, 
this model is expected to exhibit the Ostrogradski instability. The corresponding Hamilto- 
nian can be derived either by using the Dirac method |15| or directly from the Ostrogradski 
Hamiltonian as defined for higher order derivative theories. In either case, one obtains 
the expression 

2 

tt 1 \ Px 1 , 7/ 2 1 2\ 2 7 2 2 2 /„ „\ 

H PU (z,x,p z ,p x ) = — +p z x + + u 2 )x - -ui 1 oj 2 z , (2.2) 

where p z and p x are the canonical conjugate momenta of the canonical variables z and x = z 
respectively. The term linear in the conjugate momentum, p z x, represents the previously 
advertised Ostrogradski instability. At the classical level the instability may appear only 
by coupling the PU oscillator to other systems. The action fl2.1| ) leads to the fourth order 
equation of motion 

^ + K+^)g+^ = 0, (2.3) 

which is solved as superposition of two modes with frequencies ui and uj 2 in the non- 
degenerate case (c<Ji 7^ W2): 

z (t) = aie~ iu>lt + a[e iullt + a 2 e~ iul2t + a* 2 e iu " 2t . (2.4) 

From this solution we see that no instabilities in the form of exponentially growing modes 
appear in the classical solutions even though the hamiltonian has the aforementioned linear 
term associated with the Ostrogradski instability. This should not be surprising, since, 
at the classical level, only the introduction of interactions can develop the instability by 
shifting the poles of the corresponding propagator off the real axis. On the other hand, in 
the degenerate case (wi = oj 2 ) we find 

z {t) = (d + c 2 t) e~ lujt + (cj + c* 2 t) e lujt , (2.5) 

where the term linear in t comes from the fact that the roots of the characteristic equation 
for the equation of motion ( |2.3| ) are degenerate (or, in other words, the propagator has 
a double pole). Even though this term is growing in time, its growth is milder than an 
exponentially growing mode. In particular there is no imaginary propagation speed that 
would signify a classical instability. 

In the literature the PU oscillator has been already widely studied (cf. PU]). 
Here we are interested in the field version of the PU model. Let us consider the action 



5 = £ J d 4 x <p (□ + mj) (□ + mf } 



□</>) 2 - (ml + m^d^d^cp + mfm^ 2 (2.6) 



where □ = d^d^ is the d'Alembertian operator. Note that the PU scalar field <j> defined 
here is dimensionless in four spacetime dimensions. 

We now show that each spatial Fourier mode of the PU field describes a PU oscillator 
where the frequencies uj\ and lo 2 are determined by the masses mi and m 2 . From the 
action ( |2.6| ) we derive the equation of motion 

(□ + ml) (□ +m$) <p = 0. (2.7) 
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Expanding the PU field in spatial Fourier modes 



d 6 k 



(2tt) 3/2 

the equation of motion for each mode (frk becomes 



(2- 



k 4 + k 2 (m 2 + m 2 ) + m\m\ 



<Pk = 0. 



Comparing Eqs. (|2.3|) and fl2.9| ) we identify 

2 , 2 o;2 ,2,2 j 22 74, i2f 2 , 2\ , 22 

w i + w 2 = 2k + + m 2 and u^u^ = + k {m l + m 2 ) + m 1 m 2 , 
so that, as expected, the frequencies are given by: 



(2.9) 



(2.10) 



4 



-- k 2 + m\ 
k 2 +mn. 



(2.11) 



Here we see that the action ( |2.6|) describes two massive modes with masses m\ and m%. 
However, one of them is actually a ghost, sign of m 2 — m 2 determines which of the modes is 
the ghost. In the case where one of these parameters vanishes, one of the modes becomes 
massless. Notice that this corresponds to the case without the quadratic potential in the 
action, i.e., only the derivative terms remain and, nevertheless, one of the modes is massive. 
Finally, when the two masses vanish mi = m 2 = 0, we have a pure forth order theory and 
the solutions of the characteristic equation are degenerate, i.e., both modes satisfy the usual 
dispersion relation for a massless mode. This is the aforementioned degeneracy leading to 



the term linearly growing with time appearing in (2J3). From these expressions we thus see 
that each Fourier mode of a massless PU field is described by a degenerate PU oscillator 
and analogously for the massive case. 

2.2 The Stiickelberg trick 

In the following we shall show how the degenerate PU field can be identified with a 
Stiickelberg field. This field was first introduced in the Proca action for a massive vec- 
tor field to restore the U(l) gauge invariance of the theory |l(|. Thus, the vector field 
could acquire a non-vanishing mass and preserve gauge symmetry in a somewhat simpli- 
fied version of the Higgs mechanism. Here, we shall start from an action for a massless 
vector field, but with a modified kinetic term, 



S 



d 4 x 



+ -i{d^) 2 



(2.12) 



with F, 



d»A v 



d u A^ and £ and arbitrary dimensionless parameter (that can be fixed 
by the normalization of the vector field). In this action, —F 2 /4 is the usual Maxwell term, 
whereas the second term is usually introduced in the quantization of the theory as a gauge- 
fixing term and partially breaks the gauge symmetry. However, here we shall not consider 
it as a purely gauge-fixing term, but as a fully physical term. This action is no longer 
invariant under general U(l) transformations of the vector field, but only under residual 
gauge transformations An — > An + d^O which satisfy = 0. Besides the usual transverse 
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modes associated with F^ u , this action has an additional degree of freedom which appears 
here as longitudinal vector mode 1 and that is associated with d^A^. Since such a mode 
gives a scalar physical quantity, we can alternatively interpret it as a scalar mode. 



In the past [11], it has been shown that this additional mode gives a contribution to 
the cosmological energy-momentum tensor that is proportional to the metric tensor so that 
its equation of state is p = —p and it has therefore been suggested that this could play 
the role of the dark energy as an effective cosmological constant. So far, no evidences for 
longitudinal photons have been found in colliders and electromagnetism is unbroken in the 
real Universe to very high precision. Furthermore, there is no experimental evidence for 
the existence of an additional U(l) gauge field in nature. 

In this work, we show that when restoring U(l) gauge invariance using the Stiickelberg 
trick, this longitudinal mode becomes a scalar degree of freedom. Interestingly, it does not 
correspond to a normal scalar field but to the previously introduced (degenerate) Pais- 
Uhlenbeck field 



12], which has second derivatives in the action. 



Let us now restore U(l) gauge invariance by introducing a Stiickelberg field. In other 
words, we replace — > + d^cfi so that the action becomes 



S 



d 4 x 



1 

--F 



1 



1 



+ {{d^f + (D0) 2 + 2d^D0) 



(2.13) 



This action is fully gauge-invariant under gauge transformations which act on both, A^ 
and 6 



A, 



A 



n + <9 M A, 
-A. 



(2.14) 



The action (2.12) can be interpreted as the action Q2.13D in a gauge such that = 0. 
In such a gauge, the residual gauge symmetry of ( |2.12| ) remains because the condition 
\I\(p = does not fix the gauge completely, but still is invariant under gauge transformations 



6 which satisfy 06 = 0. This is analogous to the residual gauge 



symmetry that remains in standard electromagnetism after imposing the Lorenz gauge 
condition. Since we have introduced the Stiickelberg field to restore the full U(l) gauge 
symmetry, one could also consider a mass term for the vector field, as in the original 
Stiickelberg model. Then, one would obtain a massive PU field for the scalar field plus 
additional couplings to the vector field. We will not explore this possibility here because 
we are only interested in studying the massless vector field case. A way to prohibit the 
mass term (that would be allowed by the required symmetries of our action) is to impose 
the additional symmetry for the Stiickelberg field eft — > (j) + $ with # a harmonic function. 



Our action ( 2.13 ) does fulfill this symmetry, whereas the term A^d^cft that is generated 
from the mass term does not. 

The equations of motion obtained from action ( |2.13| ) are 

d v F» v + \d v A v + □</>] = , (2.15) 



□ (D<f> + dpA") = . 



(2.16) 



1 Here we refer to a 4-longitudinal mode such that A M fc M 7^ 0, as opposed with the usual terminology 
used in theories with massive vector fields where the longitudinal mode refers to the mode parallel to the 
3-momentum k. 
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Notice that the second equation is nothing but the 4-divergence of the first one. If we now 
fix to Lorenz gauge d^A^ = 0, we are left with 

□^ = ^□0, (2.17) 
n 2 <p = 0. (2.18) 

The equation for the PU field completely decouples, whereas it acts as an effective external 
conserved current term for the vector field. Such an external source is determined by the 
gradient of □</> so that it only affects the longitudinal mode of A^. In other words, the 
transverse modes completely decouple from <f> and only the longitudinal mode is affected. 
If we introduce the field = A^ — id^cj), the equations can be written as 

□£ M = 0, (2.19) 
□ 2 = (2.20) 

where now we have d^B^ = — Thus, we see that the PU field determines the longi- 
tudinal mode of the vector field B^ which, in addition, satisfies a free wave equation. In 
fact, since we still have the residual gauge symmetry, we can use it to set Bq = so that 
we obtain V • B = In this gauge, we have the free wave equations of motion for the 

transverse modes of the vector field and for □</>, whereas the longitudinal modes of the 
vector field are determined by the PU field (up to a residual gauge transformation) . Notice 



that naively imposing the Lorenz gauge condition directly in the action ( 2.13 ) leads to 



S= l -fd i x [-d^A^A" + £(D</>) 2 ] • (2.21) 



This reproduces the correct equations of motion for the transversal modes of A^ and for <fi, 
but it does not yield the correct equation for the longitudinal mode of the vector field. We 
can alternatively use a gauge such that 2 d^A^ 1 = — j^Hcf), which again leaves a residual 
gauge symmetry. Then, the equations of motion read: 

□A„ = 0, (2.22) 
\J 2 (f> = 0. (2.23) 

When imposing this gauge at the level of the action we obtain: 

(2.24) 



S= x - I <1V 



-d li A v & t A v + -^-(n<j>Y 

l + € 



which reproduces the correct equations of motion and, if we additionally impose the used 
gauge as a subsidiary condition, we also obtain the correct relation between d^A^ and 
□</>. In any case, it is more clear by looking at the equations of motion directly that 
we have completely decoupled transverse modes for the vector field satisfying free wave 
equations plus the degenerate PU field which, in addition, determines the longitudinal 
mode of the vector field. It will be useful to note that neither of the two discussed gauges 
fix the gauge freedom completely, but they leave us with the residual gauge symmetry 



2 Since the quantity d^A 1 "' + □</> is gauge invariant, we cannot impose the condition d^A^ = —□</>. For 
our choice, this only happens in the limit £ — > oo. 
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An — > + d^O, (ft — > (ft — 9 with = 0. In addition to this residual symmetry, we also 
have the aforementioned symmetry that prevents the appearance of the mass term, i.e. 
A^ — > A„, (ft — > (ft + $ with = 0. In other words, we can perform a residual gauge 
transformation with two different (harmonic) gauge parameters 9 and — 9 for A^ and (ft 
respectively. 

The transformation for the scalar field can be regarded as a generalized shift symmetry, 
since we can shift (ft not only by a constant, like in the case of a Goldstone boson, or a 
linear function of the coordinates, like in the case of the galileon field, but by an arbitrary 
harmonic function. Of course, the constant shift and the Galilean symmetry are particular 
cases of this more general symmetry. 

We know that the Maxwell term describes a well behaved theory also at the quantum 
level and with interactions. However, the lack of stability under quantization is a common 
feature in higher order derivative Lagrangians. As we have already discussed, such theories 
suffer, for instance, from the Ostrogradski instability which implies that their Hamiltonian 
is not bounded from below |m, g|. It has been argued, however, that both, the Ostrogradski 
and the ghost instability which are present in the PU oscillator can be cured and that the 
theory can be consistently quantized thanks to the presence of an unbroken VT symmetry 
that allows non-hermitian Hamiltonians to lead to unitary evolution [13]. In Section ||] we 



shall come back to this issue and develop an alternative consistent quantization procedure 
for the PU action based on the generalized shift symmetry. 

Here we also want to mention that, in the same way that we can write the Maxwell 
Lagrangian in terms of the 2-form T = dA = ^F^dx^ A dx u , as ^Maxwell = ~~ \ J J~ A *J~ 
for A = A^dx 11 , we find for the degenerate PU action, 

S PU = \i J d A x {Ucpf = 5d(ft A *5d(ft . (2.25) 

Hence d 4 x (□</>) 2 = 5d^>A*<5di^ = *d * d(ft A d * d(ft. Here 5 = — *d* is the co-differential (on 
a 4-dimensional Lorentz manifold) and *u> denotes the Hodge dual of the p-form u defined 
by 

(*w)j 1 ... in _ = —.rf^'^uji ^ , here rj is the volume form, r] = yf^gdx 1 A • • • A dx n . 

Since 55 = 0, the PU action is invariant under the transformation d(ft — > d(ft + 5T, for some 
2-form S which has the property that 5T, = d9 for some scalar field 9. The field 9 obviously 
satisfies = 5d9 = 55T, = 0. This is nothing else that the remaining gauge invariance or 
generalized shift symmetry of the PU action expressed in the language of forms. 



3. Interactions 

As we have already mentioned, it has been shown that the PU action can be quantized 
consistently. However, it happens very often that a free higher derivative theory can be 
fine but it becomes unstable once we introduce interactions. Let us do this for our theory. 
Since (ft is the Stiickelberg field of a U(l) gauge field, its natural interaction will be to 
charged fields by means of a gauge interaction. In the following we shall consider two 
explicit examples, namely a complex scalar field and a charged fermion. 
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3.1 Charged scalars 

Let us consider the action of a complex scalar field x whose interaction is mediated by our 
Stiickelberg field. We define the covariant derivative = — id^cft so that the action 
has a U(l) gauge symmetry. Under a gauge transformation x ^X an d (ft ^ (ft + 6 , the 
following action is invariant 



S 



d 4 x 



d 4 x 



D&iP^xT - v(xxl 



d^X* + XX*d^4> - id^ixd^X* ~ X*d fl X ) ~ V(XX*] 



(3.1) 



Here V is a potential which depends only on the modulus of x- If is interesting to note 
that the scalar gauge field (ft is automatically dynamical because it has derivative couplings 
to the charged field, unlike for the case of a gauge vector field. However, having a gauge 
symmetry, we can always remove it from action (|3.1[) by an appropriate gauge choice, so 
that it does not represent an actual physical degree of freedom of (|3.1[) . Nevertheless, once 
we add the free PU action to the theory, (ft can no longer be completely gauged away since 
the gauge symmetry is now reduced to gauge parameters satisfying DO = 0. Thus, the full 
interacting Lagrangian is given by 



5 



d 4 x 



d^X* + XX*d„m - id^xd^x* ~ X*d»x) ~ V( X X* 



(3.2) 



With this additional kinetic term for the gauge field (ft, it actually propagates two degrees 
of freedom, one of which can be removed using the gauge freedom so that (ft does now 
propagate one physical degree of freedom. Notice also that the addition of the PU term 
gives its corresponding propagator. 

It is also interesting to note that if the potential V(xx*) leads to symmetry breaking, we 
obtain a non-degenerate PU model from the non- vanishing vacuum expectation value (vev) 
of XX* ■ F° r instance, if we consider the usual renormalizable quartic potential V(|xl) = 
/^Ixl 2 + -^lx| 4 with fi 2 < 0, then x acquires a non-vanishing vev XX* = v 2 = — /i 2 /A and 
the quadratic term for (ft is given by 



S 



(2) 



d 4 x 



H(P(ft) 2 + v 2 d fl 



(3.3) 



Thus, we can generate the term with two derivatives analogous to the PU action with 



771,7 



777,5 



-2v /£ from a spontaneous symmetry breaking of the charged field. 



The equations of motion derived from the above action are 



DuD^x + 



dV 



dx 
dV 

dx* 



0.. 
0. 



where we have introduced the current 



f = i {x*D»x ~ XD»X*) = i (X*d»x ~ X&X*) + ^XX*d^ . 



(3.4) 
(3.5) 

(3.6) 
(3.7) 



s Here we use the fact that D M x* = (-^Vx)* so that we will not distinguish between both. 
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Notice that the first equation can be written as a conservation equation as follows: 



^(^□0-r) = O (3.8) 
so that we have a conserved charge given by 

Q = Ja 3 x (pj> - i°) . (3.9) 

Interestingly, since the x sector has a global U(l) symmetry , the current is indepen- 
dently conserved on-shell, i.e., d^j^ = so that we have the usual conservation of the 
complex field charge 

Q x = [ d 3 xj°. (3.10) 



This conservation law also implies the conservation of the current = d^Ucp that gives 
rise to the conserved charge 



Q<t> = J d 3 ;c nj> (3.11) 



associated with the PU field. The crucial fact here is that the x-dependent term disappears 
from the equation of motion of 4>. However, the equation of motion for \ still depends on 
cf>. In other words, (f> affects the dynamics of x but is itself not affected by x- It satisfies a 
"free" equation even in the presence of the coupling to x- This surprising behavior will be 
useful for the quantization of the theory because we can use the procedure of the free field 
quantization to get rid of the ghost-mode and the coupling will not reintroduce the ghost 
into the theory. 

To quantize the theory we write it in Hamiltonian form. The conjugate momenta are 
given by 

-x = f = ^V, (3-12) 

fa = 4, U 1 = ^-^ = 2 X X*<p-i(xr-XX*)-^ = f-^n<j ) , (3.13) 
o(f> at dcp 

<h = j>, n 2 = - I = ^. (3.14) 

The conjugate momentum of x* is of course 7r x * = 7r* = D°x- Thus, the Hamiltonian 
density is 

% = X^x + X^x + + ~~ ^ 

= D xD°x* - D lX D\* + V (xx*) + \^<t>? + Z ((V 2 </>) U<j> - jafy . (3.15) 

Writing this in terms of the momenta and the fields we obtain 

n 2 

w(x,^i,^2,ir x ,ni,n 2 ) = 7T X 7r*+^ 2 (x^x - xX) + + + n 2 v 2 0i 

+VxV X * + XX* (V^i) 2 - iVfa ( X Vx* - X*V X ) + V (xxl • (3.16) 
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3.2 Charged fermions 

Here we shall briefly repeat the derivations of the previous section for a coupling of & to 
a Dirac field ip. Again, by using the covariant derivative = 
charged fermion becomes 



idu<j), the action for a 



S = Jd A x [-^(-fD^ + m)^] 

d x [— ip^d^ip + iip^ipd^S — mipip] . 



(3.17) 



+ 



This action is invariant under a gauge transformation with ip e ip and (p - 
Again, having this symmetry at our disposal, the scalar field <j> can, in principle, be gauged 
away. However, when we identify it with the PU field, it actually carries two degrees of 
freedom, one of which will be physical. Thus, the full interacting Lagrangian is 



S 



d 4 x 



1 n 

-ip^d^ip + ith^d^ - rmpip + ~£ (□</>) 



This action leads to the following equations of motion 

D^ip^ — mip 
D^^ip + mip 



(3.18) 



(3.19) 
(3.20) 
(3.21) 



where we have introduced the usual U{\) current = iipj^ip, which is conserved, i.e., 
dnj^ = so that we have the usual conserved fermionic charge 



Q^ = i d 3 xipj°ip. 



(3.22) 



Again, we also have the conserved current associated with the PU field Q^. Therefore, like 
for the case of a coupling to a charged scalar field, the PU field satisfies the "free" equation, 
i.e., its dynamics is not affected by the presence of the interaction with the fermionic field, 
however the PU field <f> can, in principle, affectt the dynamics of ip. 

To end this section and for completeness, we shall compute the Hamiltonian. The 
corresponding conjugate momenta are given by 



7TV 

n x - — 

d<p 

n 2 = 



dip 

dC _ d_dC 

86 dt dip 

ins. 



Ip = 7ty 7 



iipj°ip - £D0 = j° 



(3.23) 
(3.24) 
(3.25) 



As we notice from Eq. ( 3.23 ), the field ip is proportional to the conjugate momentum -k,^ 
and so we should not consider ^ as a field like ip (cf. Jl7|]). The Hamiltonian density reads 

H = iVipifi + ni0 + n 2 <ji - C 

n 2 

= U 1( p 2 + + n 2 V 2 (/)i + 7T^7° (m + fdi) iP + (6 2 - 7°y^0i) V- (3-26) 
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In summary, in this section we have shown explicitly that, by introducing interactions 
of the PU field to charged scalars or fermions following a minimal coupling principle, 
the dynamics of the PU field is not modified. The underlying reason for this is that, 
although the quadratic term giving the free propagator for <f> only respects a residual gauge 
symmetry, the full theory still preserves the global U(l) symmetry that gives rise to the 
current conservation. Since it is precisely the divergence of the conserved current that 
enters into the equation of <fi, no effects from the charged particles on the PU fields appear. 

The above result that charged particles cannot excite the PU field can also been under- 
stood from standard electromagnetism, where fields with polarization vector proportional 
to the 4-momentum cannot be generated from conserved currents. At the quantum 
level, this is ensured by the Ward identities according to which the amplitude of any pro- 
cess involving an external longitudinal photon (in the 4-dimensional sense) vanishes. 

It is interesting to note that the above procedure to gauge a global U(l) symmetry is 
similar to making it local by introducing a longitudinal vector field. In fact, if we consider 
the covariant derivative = — iAu, nothing here imposes that the vector field 
must be transverse and only once we choose the fully gauge invariant kinetic term —F 2 /A 
for Aa, the propagating vector boson becomes transverse. This is the natural choice if we 
want the gauge boson to carry a pure massless spin-1 representation of the Lorentz group 
and also if we want to keep the full U(l) gauge invariance. However, other possibilities 
could be considered. One can, for instance, choose the kinetic term (d^A^) 2 jl for the 
vector field so that only its longitudinal component propagates. The price to pay is that 
only the residual gauge symmetry remains in the sector of the gauge boson. However, this 
is not too problematic in principle since the global U(l) symmetry remains (A^ does not 
change under a global transformation) and charge conservation is not affected. Notice also 
that, in addition to the coupling of the charged fields to 4>, we can also couple them to the 
transverse gauge field A^, in principle with a different coupling constant, i.e., they could 
be differently charged under the two different U(l) fields. 



4. Discussion 



In the previous sections we have introduced interactions for the degenerate PU field by 
following a symmetry principle, according to which all the interactions respect £7(1) gauge 
symmetry and <j> is identified with the gauge field that allows to render the gauge symmetry 
local. It is precisely this way of coupling the PU field that facilitates to consistently quantize 
the theory. We have shown that the interactions of the PU field with charged scalars or 
fermions do not modify the equation of motion for <j), although <f itself can affect the 
dynamics of the charged fields. Thus, as we shall show below, one can quantize the free 
field and use the gauge symmetry to remove the ghost-like mode. The interactions will then 
not re-introduce it. After quantizing the field, we shall discuss its cosmological relevance 
and show how it can give rise to an effective cosmological constant. 



4.1 Quantization and Stability of the PU field 



A method to quantize the PU oscillator has been proposed in [13]. This method relies on 
the fact that even if a Hamiltonian is not hermitian, it leads to a unitary quantum theory 
if it exhibits an unbroken 'P7 _ -symmetry. Here we shall not make use of this approach, but 
we shall employ the residual gauge symmetry to get rid of the undesired degree of freedom. 
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The Hamiltonian for the free PU field in terms of the conjugate momenta, 

dC _ d_dC 

06 di 86 



n. 



<06, 



dC _ 

n 2 = -tt = en? 



(4.1) 
(4.2) 



is given by 



n 2 

6 2 + -f + n 2 V 2 <; 



(4.3) 



where the expected Ostrogradski instability associated with IIi is represented by the un- 
bounded first term. This is not fatal by itself, as we have discussed above, but it can, 
in general, lead to instabilities when we couple it to another field. At the classical level, 
this instability might show up by the appearance of exponentially growing modes associ- 
ated with having an imaginary part in the propagation speed (or a propagator pole off the 
real axis). In the free theory, we have seen for the PU oscillator that actually we do not 
have such an instability, since all the modes have real propagation speeds. Only for the 
degenerate case we have a mode that grows linearly with time, associated with having a 
double pole. This represents a much milder instability than those associated with complex 
propagation speeds (imaginary parts of the propagator poles). 

On the other hand, we notice that this Hamiltonian does not have the gauge symmetry 
of the theory. However, it will obviously lead to a set of Hamilton equations that do satisfy 
the symmetry. This should not be surprising since the Hamiltonian, in general, does not 
preserve symmetries. In particular, it does not respect Lorentz symmetry, although the 
corresponding theory does. It is the set of physical observables that must respect the 
symmetries of the theory. For instance, one can always perform a canonical transformation 
that will change the form of the Hamiltonian, but leave the equations of the dynamical 
system invariant. Moreover, in the case of a field theory, it is possible to add a 3-divergence 
or a total time derivative to the Hamiltonian density without modifying the quantum 
theory [17]. In terms of the canonical variables, the gauge symmetry of the action reads 

— >• ni i2 , 6\ — > 6\ + 9, 62 — > 62 + 9 with 9 an arbitrary harmonic function. In principle, 
one could try to construct a gauge-invariant Hamiltonian by using the aforementioned 
allowed modifications, although we shall not pursue this approach here. Instead, we shall 
use the freedom given by the gauge symmetry to select a set of gauge-related modes for 
which the energy is positive. 

For our equations of motion, the field can be expanded in Fourier modes as follows: 

c\^k 1 r 

(a k + ib k kt)e i( - k - £ - kt) + (a* k - ib* k kt)e- %{k ^- 



-kt) 



(4.4) 



and a£, h\ 



1 

(2vr) 3 / 2 (2fc) 3 / 2 

After quantization, a k , b k and a* k , b* k are promoted into the operators a^, bfc 
respectively. Notice that a k is the pure gauge mode since it is modified by residual gauge 
transformations. In fact, gauge-invariant quantities (that will be determined by \D6) only 
depend on b k and b k . 

Using the above field decomposition, the Hamiltonian can be expressed as 



H 



d 3 x 



d 3 kk 



bib, 



1 



a{h k + h{ a k 



(4.5) 
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where the double dots : : denote the normal ordering operator. As anticipated, the Hamil- 
tonian is not gauge invariant which is reflected by its dependence on the gauge mode a^. 
Since we are free to perform a gauge transformation we could eliminate it from the physical 
spectrum by fixing some suitable gauge. Below, we explicitly show how to proceed. Another 
approach would be to impose an additional subsidiary condition a la Gupta-Bleuler 4 and to 
define a physical Hilbert space in which the expectation value of the gauge-dependent piece 
of the hamiltonian vanishes. The easiest way to achieve this is to define the physical states 
as those which are annihilated by the gauge mode: a^lphys) = 0. In fact, this is too restric- 
tive and can be relaxed, since imposing the most general condition (a& — iah^) |phys) = 
with a some arbitrary real number, would suffice to make the energy gauge invariant and 
positive for the physical Hilbert space. Moreover, if we only want to have positive energies, 
the parameter a of the subsidiary condition can also be complex, we just have to require 
that Im a > — 1. This can be seen by computing the expectation value of the Hamiltonian 
in a physical state: 

(H) phys = £(1 + Im q) J d 3 kk{b\b k ) phys . (4.6) 

Notice that for the energy of the physical modes to be positive for a real we need £ > and 
it has canonical energy if £ = 1. The case with Im a = — 1 in which the expectation value 
of the hamiltonian vanishes is also interesting. Then, the gauge mode exactly cancels the 
energy of the physical mode, rendering this case equivalent to the original Gupta-Bleuler 
formalism for which the scalar and longitudinal modes of the electromagnetic potential 
cancel each other . The choice of the subsidiary condition with a real seems a natural 
choice because it fully eliminates the gauge-dependent piece of the Hamiltonian and it 
amounts to identifying a physical state as the entire equivalence class of states that differ 
only by a gauge mode. In terms of the PU field, the subsidiary condition with a real 
can be written as <^ + '(io)|phys) = where is the positive frequency part of the field 
operator and to some given time. In other words, the physical states are those for which 
the expectation value of (ft vanishes at some time, i.e., (</>) p h y s = for some t = to- 
For our bi-harmonic equation of motion, one can define the scalar product 



">k,<Pk> 



i J d^x 
-i£ / d 3 x 



$fn fc ,-nf $ fe , 



(4.7) 



where T stands for transpose and we have introduced the notation 



9k 
9k 



-□0fc 
□0k 



(4.8) 



The modes that we have used to decompose the field in ( [4.4| ) are not orthonormal with 
respect to this scalar product so that the corresponding annihilation and creation operators 
do not satisfy the usual commutation relations, but they have the following commutation 
algebra: 



a/c,a 



afc,bt, 



5 {3) (k-k') 



(4.9) 



See for instance 
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with all the other commutators vanishing. The modes that diagonalize the scalar product 
are the following 5 : 



5-1/4 



>l,k 



n,k 



(2vr)3/2 (2fc)3/2 V £ 



l-y/5 



+ ikt e 



i(k-x—kt) 



1 5" 1 / 4 [2 ( 1 + \fl 



(2vr) 3 /2 (2k) 3 / 2 V £ 



+ zfci e 



i(k-x—kt) 



(4.10) 
(4.11) 



with eigenvalues +1 and —1 so that the modes 2 it have negative norm and this leads to 
a Hilbert space with indefinite metric. The PU field operator can be expanded as 



(4.12) 



A=l,2 



where now the modes are orthonormal with respect to the scalar product (4.7) and the 
annihilation and creation operators satisfy the commutation relations 



t 

\'k' 



(4.13) 



with 77 aa' =diag(l, —1). Here we see how the negative norm mode appears in the commu- 
tation relations and leads to the indefinite metric for the corresponding Hilbert space. The 
relation of these operators with at and bk is given by: 



5-1/4 



ai,fc + a 2jfc - V5(ai & - a 2 



bfc = 1/4 (ai iA . + a 2jfc 



(4.14) 
(4.15) 



Notice that the gauge-invariant mode bk is now given by (ai,fc + (i2,k)t i- e -; is given in 
terms of such a combination. The hamiltonian in terms of these operators reads: 



H 



k 



V5 L 



1 + VE) a\ fc ai )fc + ( 1 - VE) a\ k & 2 ,k + (a\ Jt a 2ifc + a!, fe ai )fc ) (4.16) 



Now we can impose the subsidiary condition (a^ — zabfc)|phys) = to obtain the physical 
Hilbert space. If we translate such a condition to the new operators it reads 



(1 — V5 — 2ia)ai jfc + (1 + \/5 - 2«a)a 2jfc 



phys 







(4.17) 



and the expectation value of the Hamiltonian in a physical state is again gauge-independent 
and positive definite. Another way of quantizing the theory that takes advantage of the ex- 
pansion in orthonormal modes is to fix the gauge such that we eliminate the negative norm 
mode and, then, quantizing the positive norm mode alone, in analogy to the quantization 
in the Coulomb gauge for standard QED. To do this, one has to fix the gauge a 2 & = a 



priori, this corresponds to choosing the gauge mode a& 



-t 1 6/ c and at 



r L bt, with 



5 Notice the interesting appearance of the golden ratio r = 1+v ^ and its inverse r 1 = — 1 ^ . 
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r the golden ratio, for the classical amplitudes. Then, the PU field operator is expanded 
in terms of the remaining positive mode, i.e.: 



d 3 k 



ai,fc<£i,fc + aj k 4>l k 



(4.18) 



Thus, these modes are orthogonal with norm +1 so that the annihilation and creation 
operators satisfy the usual commutation relations [ai^, aj . ,] = 5^(k — k'). We thus 
avoid having to work with negative norm states and negative energies from the beginning 
without having to impose a subsidiary condition. The disadvantage of this quantization 
procedure will be the lost of explicit gauge-invariance. In the quantization a la Gupta- 
Bleuler previously mentioned we work with all the modes. The expectation values of the 
physical observables are explicitly gauge-independent because they will only depend on 
h k and bt, which are the physical modes that cannot be removed by means of a gauge 
transformation. 

The natural concern arising with the discussed procedure to identify the physical states 
is whether the introduction of interactions will spoil it. However, as we have shown in the 
previous section, if interactions are introduced following a symmetry principle, this will 
not be the case because, while the PU field will affect the sector of charged particles, 
its own equation of motion remains unaffected. The reason for this is that introducing 
the couplings as U(l) gauge interactions leads to the presence of the associated Noether 
current, whose divergence is precisely the new term appearing in the equation of motion 
of the PU field. Thus, the conservation of this current also implies that the equation of 
the PU field is not modified. This is the crucial point to guarantee that the discussed 
procedure to identify the physical states remains unaffected when introducing interactions 
because the field decomposition given in ( |4.4| ) is also valid (and exact) in the presence of 
interactions with charged particles. 

4.2 Cosmological relevance 

Let us finally show how the PU field can actually play the role of dark energy. For that, 
we first note that □</> satisfies the usual equation for a massless scalar field which, for a 
homogeneous field in a FLRW universe with metric ds 2 = dt 2 — a(t) 2 dx 2 , reads 



dT 2+m dt 



□0 = 0. (4.19) 



The solution of this equation is given by 

□0(t) =C X + C 2 J^ (4.20) 

where Ci^ are integration constants. Since the CVmode decays throughout the expansion 
of the universe, only the constant mode C\ is relevant at late time. This is simply the 
well known result that a massless 6 scalar field is frozen on super-Hubble scales. If we now 
compute the energy-momentum tensor for the PU field we obtain 

= ig^ + dx^n^j - 2Hd { ^d u p<p. (4.21) 



'More precisely, a scalar field with a mass much smaller than the Hubble expansion rate. 
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Thus, it becomes apparent why the PU field can drive an accelerated expansion. At late 
times, □</> is constant so that we obtain that the above energy-momentum tensor is simply 



V = i(D0) 2 ^. (4.22) 

This is the form of the energy-momentum tensor of a cosmological constant with the value 
A = 8-7rG^(D0) 2 . Notice that we must have £ > for the effective cosmological constant to 
be positive. This is also what one would expect since, if □</> is constant, it contributes as 



a cosmological constant in the action. This result is in agreement with Ref. [19| where the 
cosmology of a theory whose action is a general function of the D'Alembertian of a scalar 
field is explored and it is found that de Sitter is an attractor for this type of theories. This 
is also expected from the results of Refs. [11], where the symmetry breaking term (V^A 11 ) 2 
gives rise to the appearance of an effective cosmological constant on super-Hubble scales. 
This relates to our result here because our PU field plays the role of the additional degree 
of freedom introduced by the symmetry breaking term there. 

Since the field □</> behaves like a massless scalar field, the primordial power spectrum 
generated from its quantum fluctuations during a de Sitter inflationary phase will be the 
usual scale invariant one given by 7 Pas — Hj where Hj is the constant Hubble parameter 
during inflation. Of course, in a more realistic quasi de Sitter inflationary phase, we expect a 
slightly tilted power spectrum with spectral index proportional to the slow roll parameters. 
However, the de Sitter expression will suffice for our purpose. 

Interestingly, the scale of the effective cosmological constant is set by the scale of 
inflation M 2 A ~ (D^) 2 ~ Hj ~ (M|/M p ) 4 , with Mi and M p the inflationary and Planck 
scales respectively and we have used the Friedmann equation to relate Hi with Mi. Since 
we know that A ~ Hq where Hq denotes the Hubble parameter today, we obtain Mj ~ 
HqMj } . If we use the corresponding values for Hq and M p we find Mj ~ ITeV, i.e., 



the electroweak scale. This is the result that was also found in Ref. [11] where the role 
of effective cosmological constant is played by (V^A' 1 ) 2 . Also, in [20] it is shown that 



quantum fluctuations of a scalar field during inflation can produce dark energy provided 
its mass is smaller than the Hubble expansion rate today. 

This is a general feature of this type of models: If we have an action without any 
dimensionfull parameter that effectively gives a massless scalar field and whose energy 
density is constant on super-Hubble scales, then the value of the effective cosmological 
constant that is generated during an inflationary phase taking place at the electroweak 
scale coincides with the observed value. However, the scalar field must arise from some 
non-standard mechanism, since the energy density of a standard scalar field is diluted 
by the expansion of the universe. In the present case, the effective scalar field is the 
physical degree of freedom remaining in the degenerate PU model that is associated with 
the d'Alembertian of the field, while in [|ll]] its role is played by the divergence of the vector 
field. 



7 The power of Hf appearing here as opposed to the power Hj found for the usual scalar field can be 
understood from dimensional arguments, but the underlying reason is that (f> satisfies a fourth order equation 
so that it has a different normalization. Thus, although \3<j> satisfies the same equation as a massless scalar 
field, the power spectrum amplitude is different. 
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5. Conclusions and outlook 

In this work we have shown that the degenerate Pais-Uhlenbeck field arises naturally as 
the Stiickelberg field that restores the full U(l) gauge symmetry in the action for a vector 
field including the usual Maxwell term plus a term oc {d^A^) 2 . We have seen that, in a 
given gauge, the Stiickelberg field decouples from the vector field and satisfies the equation 
of motion of the degenerate PU field. Moreover, it determines the value of the longitudinal 
mode of the vector field so that it can be identified with it. After fixing to the decoupled 
gauge, we are left with the residual gauge symmetry analogous to the one remaining after 
imposing the Lorenz condition in standard electromagnetism. For the scalar field, this 
residual symmetry is an invariance under the addition of an arbitrary harmonic function 
so that it is a generalized version of the shift or galilean symmetries. 

Although we have a fourth order theory, it propagates only one physical degree of 
freedom, since the second degree of freedom can be removed by the residual symmetry. 
Also the ghost-like mode that one expects from a fourth order theory can be removed 
from the physical spectrum by an appropriate gauge choice so that the free Hamiltonian 
only contains the gauge-invariant and positive energy of the physical degree of freedom. 
Although this can easily be achieved in the case of the free field, the introduction of 
interactions could spoil the mechanism by exciting the ghost-like mode. However, we have 
shown that if the couplings are introduced in a manner that the PU field acts as the gauge 
boson associated with a U(l) symmetry, then the excising mechanism is not spoilt. 

We have explicitly worked out two particular cases: coupling to a complex scalar field 
and to a Dirac fermion. We have seen that in both cases, the equation of motion for the PU 
field remains unaffected thanks to the current conservation granted by the U(l) symmetry. 
This is crucial for the stability of the theory because guarantees that the ghost will not 
be reintroduced in the theory by the interactions and the quantization procedure remains 
consistent even in the presence of interactions. 

Even though particles that are charged under such a U{\) group are affected by the 
presence of the PU field, they cannot generate the PU field. The only way in which it can 
be produced is from quantum fluctuations during the inflationary era in the early universe. 
Once the quantum fluctuations are amplified, the super-Hubble modes contribute as an 
effective cosmological constant to the energy-momentum tensor whose scale is determined 
by the scale of inflation. The homogeneous evolution is exactly the same as in standard 
ACDM. However, as the cosmological term is truly the dynamical PU field, it can be 
perturbed and the evolution of the cosmological perturbations including it will differ from 
standard ACDM. This provides us with a mechanism that can help to discriminate this 
model from a pure cosmological constant. 
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